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Abstract. Sunspots, faculae and the magnetic network contribute to solar irradiance variations. The contribution due to faculae 
^ ' and the network is of basic importance, but suffers from considerable uncertainty. We determine the contrasts of active region 

faculae and the network, both as a function of heliocentric angle and magnetogram signal. To achieve this, we analyze near- 
simultaneous full disk images of photospheric continuum intensity and line-of-sight magnetic field provided by the Michelson 
Doppler Interferometer (MDI) on board the SOHO spacecraft. Starting from the surface distribution of the solar magnetic field 
we first construct a mask, which is then used to determine the brightness of magnetic features, and the relatively field-free part 
of the photosphere separately. By sorting the magnetogram signal into different bins we are able to distinguish between the 
contrasts of different concentrations of magnetic field. We find that the center-to-limb variation (CLV) of the contrast changes 
strongly with magnetogram signal. Thus, the contrasts of active region faculae (large magnetogram signal) and the network 
■ (small signal) exhibit a very different CLV, showing that the populations of magnetic flux tubes that underly the two kinds of 

features are different. The results are compatible with, on average, larger flux tubes in faculae than in the network. This implies 
that these elements need to be treated separately when reconstructing variations of the total solar irradiance with high precision. 
We have obtained an analytical expression for the contrast of photospheric magnetic features as a function of both position on 
"qJj , the disk and spatially averaged magnetic field strength, by performing a 2-dimensional fit to the observations. We also provide 

a linear relationship between magnetogram signal and the fi = cos ((9), where 9 is the heliocentric angle, at which the contrast 
is maximal. Finally, we show that the maximum contrast per unit magnetic flux decreases rapidly with increasing magnetogram 
signal, supporting earlier evidence that the intrinsic contrast of magnetic flux tubes in the network is higher. 

><. 
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1. Introduction 

Radiometers on board satellites launched during the last three 
decades (NIMBUS-7, SMM, UARS, EURECA, SOHO) have 
revealed that the total solar irradiance, also referred to as the 
solar constant, changes on a variety of time-scales. Solar irra- 
diance variations on scales of days up to the solar activity cycle 
length are closely related to the evolution of the solar surface 
magnetic field, because the emergence and evolution of active 
regions (AR) on the solar surface is reflected in the irradiance 
records (Lean et al. |1998| Fligge & Solanki |2000b[ ), Sunspots 
and active region faculae are considered to be the dominant 
contributors to solar irradiance changes on time-scales of days 
to weeks. Space-based irradiance records have also established 
a variation of about 0.1% of the irradiance in phase with the 
1 1 year solar activity cycle, giving as a result a brighter Sun 



around activity maximum (Chapman 1987; Willson & Hudson 
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1988). The origin of the long-term increase of the irradiance 



between activity minimum and maximum is still widely de- 
bated. It is expected that small-scale magnetic elements that 
compose the enhanced and quiet network contribute substan- 
tially to the observed irradiance increase during activity max- 



imum (Foukal & Lean 1988; Solanki & Fligge 2001; Fligge 



& Solanki 2000a). Nevertheless, other mechanisms of non- 



magnetic origin have also been proposed, based, for example, 
on temporal changes in the latitude-dependent surface temper- 



ature of the Sun (Kuhn et al. 1988). Other authors have tried to 



explain these variations by modelling structural change s in th e 
convection zone during the solar cycle (Balmforth et al. 1996) 



The photospheric magnetic field is concentrated in discrete 
elements whose diameters range from less than a hundred to 
several tens of thousands of kilometers. The brightness signa- 
ture of these magnetic features is a strong function of their he- 
liocentric angle and their size; sunspots are dark while small 
flux tubes are ge nerally brig ht; faculae appear brighter near 
the limb (Solanki 1993 2001). However, our knowledge of the 
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brightness of small scale magnetic features, groups of which 
form faculae and the network, is incomplete (e.g. Solanki 



1994) 



Faculae, bright structures seen in the photosphere cospa- 
tially with chromospheric plages, are associated with magnetic 
fields. At high resolution, they consist of many unresolved 
small continuum bright points, with diameters of about 100 km, 



called facular points (Muller 1983; Bergeret al. 1995). The ob 



served zoo of small magnetic features is unifyingly described 
by the concept of the small flux tube. To predict their physi- 
cal properties, different models for small flux t ubes have been 
constructed (e. g., Sp ruit |l97rj| ; Deinzer et al. |l984aj |l984b ; 



Knolker et al. 1988; Knolker & Schussler 1988; Grossmann 



Doerth et al. |1989| ; Steiner et al. |1996| ). According to this model 
faculae are conglomerates of evacuated flux tubes with hot 
walls and a hot or cool floor (corresponding to an optical depth 
of r = 1) depending on the evacuation and diameter of the flux 
tube. The model predicts a certain CLV of the contrast for a 
particular diameter of the underlying flux tubes. 

This model assumes that inside each small flux tube the 
magnetic field is of the order of a kilogauss, but practically 
zero outside. Due to the magnetic pressure the flux-tube inte- 
rior is evacuated, so that r = 1 is reached along its walls (which 
are bright due to radiation leaking in from the surroundings, 
i.e., there is a horizontal flux of energy into the tube). Flux 
tubes can be dark at disk center if suppression of convective 
energy transport within the tube is included (e.g., Spruit 1976; 
Deinzer et al. |1984a[ |1984b| ; Knolker et al. |1988[ ), resulting in 
a cooling of the deeper layers. When the tube is sufficiently 
broad, the horizontal optical depth between the wall and the 
tube center is large and most of the radiation cannot reach the 
center. In this case, the tube floor remains dark at its center. 
But if the tube is sufficiently slender, the horizontal flux of en- 
ergy can reach the center of the tube; then, the interior of the 
tube is heated, the vertical energy flux increases, and the tube 
turns bright, even at disk center. In this scenario, the transi- 
tion between sma ller br ight points and larger dark micropores 
(e.g., Topka et al. 1997 ) occur s at a diameter of about 300 km 
(e.g., Grossmann-Doerth et al. 1994 ), and therefore micropores 
would fill the gap between small bright points and larger dark 
pores. Micropores are predominantly found in active regions, 
while bright points are the main constituents of the network. 
When observed near the limb, the heated walls of the tube be- 
come visible, and therefore the contrast increases. 

Observations also provide evidence that the contrast, as 
well as the underlying thermal structure, depends on the size 
of the flux tubes (e.g., Keller 1992), but more commonly on 
the strength of the magnetogram signal (Frazier 1971; Spruit 
& Zw aan|l98~ll ; Solanki & Stenflo |983 Solanki J§8§ Zayer 
et al. |1990|; Solanki & Brigljevic |1992fc Topka et al. |1992, 



1997; Lawrence et al. 1993 Grossmann-Doerth et al. 1994), 



and hence a test of flux tube models is possible with such data. 
However, since most flux tubes are not resolved, it is necessary 
to have a well-defined and constant spatial resolution of the 
observations with which to compare the models. With ground- 
based data, the basis of practically all facular contrast CLV 
measurements to date, this criterion is hard to meet. 



Differences in spatial resolution, caused by variable see- 
ing, may indeed partly explain the variety of measured con- 
trast CLVs (e.g., Libbrecht & Kuhn |l984[ Unruh et al. |l999|) 



Libbrecht & Kuhn |1984fc Unruh el 
Other possible factors are differences in wavelength, spectral 
resolution and the magnetic filling factor of the observed fea- 
tures (Solanki 1994 ). The problem posed by variable seeing can 
be circumvented by employing data recorded in space, while 
the magnetic filling factor can be estimated with the help of 
cospatial and cotemporal magnetograms. Only relatively few 
contrast investigations including the magnetogram signal can 
be fo und in the lit erature (e.g. Frazier 1971 ; F oukal & Fowler 



1984; Topka et al. 1992, 1997; Lawrence et al. 1993). 



Here we add to this list using data fr om the MDI instrument 
on board SOHO (Domingo et al. 1995 ); their main advantages 
are: 

— seeing effects due to the Earth's atmosphere are avoided, 

— measurements are made in the pure continuum (i.e. spectral 
resolution is not a problem), 

— the 20-minute averaged MDI magnetograms have a reason- 
ably low noise level, 

— a large and homogeneous data set is available, 

— the characteristics of the instrument and the data sets are 
well known and stable, 

— magnetograms and intensity images are obtained regularly 
by the same instrument with exactly the same spatial reso- 
lution, so that a one-to-one identification of brightness with 
magnetogram signals can be made. 

The main disadvantages of the MDI data are: 

— a relatively low spatial resolution with a pixel size of 2 x 2", 

— measurements are available at only a single wavelength. 

The purpose of this paper is to present new high-quality 
measurements of the contrast of the photospheric bright fea- 
tures as a function of both heliocentric angle and magnetogram 
signal and to obtain an analytical function that predicts their 
contrast given a position on the disk and a magnetic signal 
value. Such measurements are expected to be of use not only 
to constrain models of flux tubes, but also to im prove the mod- 
elling of the solar irradiance (Lean et al. 



1998). Uncertainties 



in the contrast of faculae and the network are one of the ma- 
jor sources of error in the modelling of solar irradiance varia- 
tions. Employing MDI data to obtain the contrast as input for 
irradiance modelling is of particular interest since MDI magne- 
tograms have already been succesfully used for such modelling 
(Fligge& Solanki |2000a| ). 

In Sect. H we present the data sets used and the analysis pro- 
cedures. In Sect. |we describe the results, which are discussed 
in Sect. [|. Finally, our conclusions and a summary are given in 
Sect. g. 

2. Data and analysis procedure 

2.1. Data sets 

The Solar Oscillations Investigation/Michelson Doppler 
Imager (SOI/MDI) instrument is a state-of-the-art helioseis- 
mology experiment and magnetograph on board the SOHO 
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spacecraft, devoted to study the interior structure and dy- 
namics as well as the surface magnetic field of the Sun. This 
instrument gives an image of the Sun on a 1024 x 1024 CCD 
camera, and can observe in two spatial resolution modes, full 
disk and high-resolution of the central part of the disk (HR). 
We are interested in the full disk measurements, which have a 
field of view of 34 x 34' and a pixel size of 2 x 2". Two tunable 
Michelson interferometers allow MDI to record filtergrams 
centered at five wavelengths across the Ni I 6768 A absorption 
line. From the filtergrams, MDI computes the following six 
observables: Doppler velocity, continuum intensity, line depth, 
longitudinal magnetic field, horizontal velocity and limb 
position. The SOI/M DI instrument is described in detail by 
Scherrer et al. ( |l995| ). 

The products of interest for our work are the full disk mag- 
netograms and continuum intensity images. Magnetograms 
only measure net magnetic flux per resolution element, there- 
fore the signal is not the true magnetic field strength B, inside 
a flux tube, but its longitudinal component, (|B| cos 7}, aver- 
aged over the pixel, where 7 is the angle between the mag- 
netic vector and the line of sight. For simplicity, we hereafter 
refer to (|B| cos 7} as B. In a 2-component model of the mag- 
netic field, with magnetic flux tubes of field strength B cov- 
ering a fraction a of the solar surface separated by a field- 
free component covering (1 — a) of the surface, we can write 
(|B| C0S7) as a|B| COS7. Since the true field strength |B| lies 
in a relatively narrow range of 1000-1500 G for all magnetic 
features except intranetwork elements (Solanki et al. 1999), 
and cos 7 w cos w // is a reasonable approximation (0 is 
the heliocentric angle), the strength of the magnetogram signal 
mainly provides information on the magnetic filling factor a. 

MDI magnetograms are usually obtained every 96 minutes, 
with the exception of periodic campaigns in which 1 -minute 
cadence measurements are available. The 1-er noise level for a 
one-minute longitudinal magnetogram is 20 G. Full disk con- 
tinuum intensity images are taken each minute with a noise 
level of 0.3%. 

The analyzed data set consists of nearly simultaneous mag- 
netograms and continuum intensity images recorded on 10 days 
in the period February to October, 1999, as shown in Table [j]. 
The time of the observations is given for the averaged magne- 
tograms (see Sect. 2.2). It corresponds to the middle of the 20 
minute integration time. These days were chosen because they 
belong to a high activity period so that everything from quiet 
network to intense active-region plage was present on the solar 
surface. The sample contains active regions spread over almost 
all fj, — cos 8 values. They are also generally well separated in 
time in order to avoid duplication. 



2.2. Reduction method and analysis 

We employ averages over 20 single magnetograms, taken at a 
cadence of 1 per minute, in order to reduce the noise level suf- 
ficiently to reliably identify the quiet network. The individual 
magnetograms were rotated to compensate for the time differ- 
ence before averaging. Intensities are standard 1 -minute im- 
ages. Care has been taken to use intensity images obtained as 



Table 1. Selected days and times (hours/min/sec) during 1999 
at which the averaged magnetograms analyzed here were 
recorded (see text for details). 



1999 observation dates Time (UT) 


February 13 


00:10:02 


February 20 


04:10:02 


May 14 


00:10:03 


May 28 


06:10:03 


June 25 


01:10:02 


July 2 


03:10:02 


July 10 


01:10:02 


August 7 


00:10:02 


October 12 


09:10:03 


October 15 


06:10:03 



close in time to the magnetograms as possible. In all cases but 
one, the two types of images were recorded within 30 min- 
utes of each other, with 37 minutes being the highest differ- 
ence. The intensity images have been rotated to co-align them 
with the corresponding average magnetogram. Intensity images 
have also been corrected for limb-darkening effec ts usin g a fifth 
order polynomial in // following Neckel & Labs ( 1994 ). Our fi- 
nal data sets are pairs of co-aligned averaged magnetograms 
and photospheric continuum intensity images for each of the 
10 selected days. Both types of images can be compared pixel 
by pixel. An example magnetogram and the corresponding in- 
tensity image recorded on October 12, 1999 are shown in Fig. [j] 
(top and middle panels). 

We have determined the noise level of the MDI magne- 
tograms and continuum images as a function of position over 
the CCD array. The standard deviation for the magnetic sig- 
nal has been calculated using a running 100 x 100 pixel box 
over the solar disk, with the exception of the limbs, which 
were avoided by masking out an outer ring of 75 pixels width. 
This process was applied to several 1996 low activity magne- 
tograms, in order to avoid artifacts introduced by the presence 
of active regions. After that, their median was determined to 
eliminate the possible remaining activity. A second order sur- 
face was then fitted to the result and extrapolated to cover the 
whole solar disk. The resulting noise level, <7 mag , shows an in- 
crease towards the SW limb that probably includes some ve- 
locity signal leakage. In Fig. || we show the calculated stan- 
dard deviation for the 20-minute averaged magnetograms. Note 
that when applying this noise surface to our data we have as- 
sumed that the MDI noise level has remained unchanged be- 
tween 1996 and 1999. 

A similar procedure has been used to determine the mean 
and standard deviation of the quiet Sun continuum intensity 
for each selected day, (7 qs ) and cri qs respectively, where the 
subscript qs denotes "quiet Sun". Every pixel in the running 
mean box with an absolute magnetic signal value below 0.5 
times cr mag (i.e. pixels with corresponding magnetogram signal 
between approximately —2.5 and 2.5 G) has been considered 
as a quiet Sun pixel. 

The surface distribution of solar magnetic features that pro- 
duce a bright contribution to irradiance variations, is identified 
by setting two thresholds to every magnetogram-intensity im- 
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Fig. 2. Standard deviation (in G) of the 20-minute averaged 
MDI magnetograms. An increase of the noise in the direction of 
the SW limb is evident. The shadowed contours indicate some 
of the values. Both the surface and the contours represent the 
standard deviation. 



■ ■■■■ ■'■ ' :. "'. 



Fig. 1. Example of a 20 minute averaged MDI magnetogram 
(top panel), the corresponding intensity image after removal of 
limb-darkening (middle panel) and the resulting contrast mask 
(lower panel) for October 12, 1999. 



age pair. The first threshold looks for magnetic activity of any 
kind, and is set to ±3cr mag , which corresponds, on average, to 
15 G. As we are only interested in bright magnetic features, 
the second threshold masks out sunspots and pores by setting 
all pixels with a continuum intensity 3cri qs below the average 



to a null value. To reduce false detections, even at the risk of 
missing active pixels, we reject all isolated pixels above the 
given thresholds assuming that they are noise. 3 10 4 out of 10 7 
analyzed data points are rejected in this way. After this step, 
we find that 6% of the pixels satisfy both criteria. Using both 
thresholds we construct a mask of the contrast of bright features 
for each day. The result of applying the mask derived from the 
magnetogram (top panel) and intensity image (middle panel) 
shown in Fig. [j], is displayed in the bottom panel of that fig- 
ure. Note that only features that lie above the given thresholds 
in the magnetogram and the intensity image are indicated by 
white pixels. Sunspots near the NE limb, for example, do not 
appear in the mask, but faculae surrounding those sunspots are 
well identified. Smaller features belonging to the magnetic net- 
work are also pinpointed outside of the active regions, although 
weaker elements of the network may well be missed. For each 
pixel with coordinates (x,y), the contrast Cf ac is defined as: 



Cf ac (x,y) = 



I(x,y) - (I qs )(x,y) 
(I qs )(x 7 y) 



(1) 



Contrast, magnetic field strength averaged over the pixel 
and position, represented by the heliocentric angle fj, = cos(0), 
are calculated for each selected pixel. Finally, for each of these 
parameters the pixels above the thresholds for each of the 10 
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selected days are put together into vectors of about 6 10 5 ele- 
ments, which should provide adequate statistics for a detailed 
study of the facular and network contrast. 

The method used in this work resembles that employed by 
Topka et al. (1992 1997), although our magnetic threshold is 
much lower due to the less noisy magnetograms used. The an- 
gular resolution, however, is also considerably lower, but it is 
constant. 

3. Results 

We have analyzed the AR faculae and network contrast depen- 
dence on both /i and the measured magnetic signal, B. It is im- 
portant to recall that the observed magnetic signal drops to zero 
at the limb, even if strong magnetic field regions are present. 
This is a straightforward consequence of the fact that magne- 
tograms are only sensitive to the line-of-sight component of the 
magnetic field and that the magnetic field is mainly vertical. To 
compensate this effect to first order we have worked with B/fi 
(i.e. (|B| cos7)/cos60 instead of B. Second order effects due 
to radiative-transfer effects or finite thickness of flux tubes re- 
main. 

We have binned the B//i values into eight intervals that 
range from the threshold level of, on average, 15 G to about 600 
G. The intervals have been chosen so that they roughly contain 
the same number of pixels each, with the exception of the high- 
est Blii bins which have fewer pixels. By sorting the magnetic 
field strength into different bins we can distinguish between the 
CLV of the magnetic features present in regions with different 
filling factor a. Fig. || displays the contrast, Cf ac , as a function 
of /i, individually for every B//i interval. A second order poly- 
nomial has been fitted (employing least squares) to guide the 
eye and a dashed line indicating Cf ac = has been included 
for clarity. To avoid overcrowding we have binned data points 
in sets of 40 before plotting. 

Fig. H reveals a clear evolution of the behaviour of the con- 
trast from one B//i interval to another. Network features (top 
left panel) show a low and almost constant contrast, as com- 
pared with the very pronounced CLV of the contrast for active 
region faculae (bottom panels). Intermediate cases show a pro- 
gressive increase of the contrast towards the limb as well as an 
increasingly pronounced CLV. When B/fi< 200 G, the con- 
trast peaks at ii ~ 0.5, while for magnetic signals Blfi > 200 
G this maximum shifts to lower values of [i (see Fig. || and its 
discussion in the text for a more quantitative analysis). Note 
also that for Blfi > 200 G the contrast is negative around disk 
center, while it is positive in the network (i.e. for smaller Bj [£). 
We will return to this point in Sect. The large fluctuations of 
the contrast near the limb for intermediate and high magnetic 
signals are due to the distribution of active regions on the ten 
selected days. 

Fig. |] shows the contrast as a function of B/fi, for different 
positions on the solar disk. The solar disk has been divided into 
eight bins of /.i, centred on ji values ranging from 0.96 (disk 
center, top left) to 0.3 (limb, bottom right). Note that to keep 
the number of points in each bin approximately equal, /i-bins 
lying closer to the limb are wider than the ones around /i = 
1, showing therefore some overlap. Each point of the figure 



is obtained by binning together 40 points of data with similar 
B//i, and second degree polynomials have also been plotted, as 
in Fig. |. 

This figure shows that, in general, the contrast initially in- 
creases with Blii before decreasing again. We expect it to con- 
tinue decreasing for even larger values of B/fj, representing 
pores and sunspots. At large fi the initial increase is small and 
the contrast basically decreases with B/fi, while at small /i it 
mainly increases. For fi = 1, points with large B/fi show a neg- 
ative contrast (as in the lower panels of Fig. ^J), while points at 
the limb always have positive contrasts. The bins with fi > 0.82 
do not display data for high magnetic signals, because their in- 
tensity is below the intensity threshold. 

Given the regular behaviour of the contrast as a function of 
[i and magnetogram signal, Cf ac (/z, Bj fi), it seems appropiate 
to search for an analytical expression for this dependence. We 
have performed a multivariate analysis using a (/i, B/fi) grid. 
The fi values have been binned linearly, with A/i = 0.1. B/fi 
bins have been chosen to be equally spaced on a logarithmic 
scale, with Alog(_B//i) = 0.05, in order to compensate for 
the fact that magnetic signals are mostly concentrated towards 
lower values (Figs. || and |]). The dimensions of the grid are 
0.1 < fi < 1 and 17 G < (B/fi) < 630 G, resulting in a 9 x 
31 bins grid. We do not consider points with Blfi > 600 G 
to exclude bright features that might belong to pores observed 
near the limb. We are aware of the fact that 600 G is an arbitrary 
value for such a cutoff. 

Each bin of the grid is defined by the averaged values of 
the contrast, fi and Blfi, over all the data points of that bin. 
Although the curves in Fig. ||are only intended to guide the eye, 
they do reveal that second order polynomials fit the contrast as 
a function of fi well. We have fitted the bidimensional array 
of contrasts by a second order polynomial function of fi and a 
cubic function of Blfi of the form: 



c^B/^^Yt^H-y, 



(2) 



where i runs from 1 to 3, j runs from to 2 and djj are the 
coefficients of the fit. The result of the fit is a surface of sec- 
ond order in position fi and third order in magnetic signal. The 
coefficients of the multivariate fit a, ; are: 



C iac {fi,B/fi) = 10~ 4 [0.48 + 9.12/1-8. 50/i 2 ] 



(3) 



B 
fi 

10~ 6 [0.06 - 2.00fi + 1.23m 2 ] (j; 
1Q- 10 [0.63 + 3.90^ + 2.82^] 



The terms of Eq. (|2|) are grouped in Eq. (|3]) to make clearer 
the quadratic dependence of Cf ac (/i, B/fi) on fi and the cubic 
dependence on B/fi. When Blfi is small (< 100 G), the first 
order term in Blfi provides the dominant contribution. When 
Blfi is large (> 200 G), first and second order terms in B/fi 
dominate the contrast, modulated by the contribution of the cu- 
bic term which plays a role in this range of magnetic signals. At 
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Fig. 5. Polynomial surface of second order in fi and third order 
in B/fi obtained from a multivariate fit performed to the grid 
of contrasts, covering fi and B/fi values. Dashed vertical lines 
project the corners of the plotted surface onto the fi-B/fi plane 
and indicate the region spanned by the fit. 

disk center (/i = 1), those terms result in a dominant negative 
contribution (Fig. g). Note that the contrast is constrained to go 
through zero when B/fi=0, as expected for the quiet Sun. 

The best-fit surface is shown in Fig. ||J The grid corresponds 
to the linear /i-bins and the logarithmic B/fi-bins taken for the 
fit. The shape of this surface is quite congruent with that shown 
by the observed contrast in Figs. || and [I| Note that the function 
given in Eq. (^|) is valid only for the wavelength and spatial 
resolution of the MDI data, 6768 A and 2" , respectively. For 
other values of these parameters we expect another dependence 
on fi and B/fi. In particular, the absolute value of the contrast 
is expected to change. 

To better estimate how this analytical surface fits the be- 
haviour of the measured contrasts, we have sliced the surface in 
both directions, fi and magnetograph signal, and then compared 
the result with the measured values. In Fig. || the fitted surface 
is sliced along the /i-axis (solid curve), at three sample mag- 
netic signal ranges, representative of low (top panel), medium 
(middle panel) and high (lower panel) B/fi values. Dots rep- 
resent measured contrasts. To avoid very crowded plots each 
dot represents 250 (top panel), 100 (middle) and 25 (bottom) 
data points, respectively. The different amounts of binning re- 
flect the non-uniform distribution of points over the B/fi range. 
The multivariate regression surface fits quite well the plotted 
dependence of the contrast, although minor deviations are vis- 
ible at small B/fi. Fig. shows slices of the modeled surface 
along the B/fi axis (solid curves) and the corresponding binned 
data (dots), at three sample positions on the solar disk, from 



disk center (top panel) to near the limb (lower panel). The fit- 
ted curves now deviate somewhat more from the data points, 
most significantly for 0.54 < fi < 0.66, where the discrepancy 
can reach 0.01 in contrast. 

The multivariate analysis yields to an expression for the 
contrast of photospheric bright features, Cf ac (/U, B/fi), that 
cannot be directly compared with previous studies because, 
to our knowledge, no similar work has been done before. 
Quadratic f unction s have already been used by other authors 
(e.g. Foukal 1981 ) to fit the dependence of the facular contrast 
on position over the disk, although most of them use a func- 
tion of the form Cf ac (fi) — b(l/fi — a) (Chapman 1980). A 
quadratic function agrees quite well with the CLV proposed by 
the hot wall model. A cubic function has been used for fitting 
the dependence of the contrast on magnetic strength. In this 
case we do not have a physical reason, only the goodness of the 
fit with respect to other bivariate functional dependences tried 
(see Figs. || and ^) and the requirement to force the contrast 
through zero for a disappearing magnetic signal. We suspect 
that, in order to obtain a better empirical description of the de- 
pendence of facular contrast on fi and B/fi, a larger number of 
free parameters is required. 

The dependence of the peak of Cf ac on B/fi is shown 
in Fig. H The /i-values at which Cf ac peaks, ^ max , have 
been represented against the corresponding magnetic signal in 
Fig. pk. Fig. ||b shows the peak Cf ac values reached, Cf™ x (see 
Fig.Bh, plotted as a function of B/fi. Finally Fig. Kp shows 
C"™ x I (B I fi) plotted against B/fi or, in other words, the de- 
pendence on the magnetic signal of the specific contrast per 
unit of magnetic flux. Errors in fi are estimated from the differ- 
ence between the peak of the best-fit curves in Fig. || and the 
peak obtained directly from the data points. Error bars in B/fi 
correspond to the size of the B/yU-intervals used in Fig. ||. 

A linear regression adequately describes the dependence of 
Mmax on B / fi for the precision achievable with the current data. 
The best fit straight line (solid line in Fig. pk) is 



5.60 10" 1 -4.97 lQ- 4 (B/fi). 



(4) 



Thus, /Umax = 0.56 ± 0.02 when B/fi tends to zero, and 
Mmax = for B/fi w 1120 G from extrapolations of this curve. 
Fig. ^c implies that the contrast per unit of magnetic signal de- 
creases strongly with increasing magnetogram signal. Since in- 
dividual flux tubes are not resolved by MDI, we cannot infer the 
intrinsic contrast of a flux tube from Fig. ||b, which obviously 
shows the same pattern as Fig. [|. However, by normalizing by 
B / fi we obtain a quantity that is roughly proportional to the 
intrinsic brightness of the flux tubes (assuming that the field 
strength of the elemental magnetic flux tubes lies in a narrow 
range as mentioned in Sect. ^). 

4. Discussion 

4. 1. Comparison with previous observations 

Comparison with other contrast observations is not easy be- 
cause of the differences in the selected wavelength, spatial res- 
olution, range of studied heliocentric angles, magnetic filling 



A. Ortiz et al.: On the contrast of solar photospheric faculae and network elements 



7 



50S(B/ / u)<80 



200S(B/|i)<300 




400S(B/»<500 




0.0 



0.2 



0.4 0.6 



0.8 



Fig. 6. Comparison of cuts through the surface (solid curves) 
returned by the multivariate analysis and the measured con- 
trasts (dots) as a function of /i, for 3 sample bins of corrected 
magnetic signal. Every dot represents 250 (top), 100 (middle) 
and 25 (bottom) data points. 
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Fig. 7. The same as Fig. |(| but for cuts along the B/fi-axis (solid 
curves) made at three positions on the solar disk. Dots represent 
measured contrasts. The plotted curves represent the same fj, 
ranges as those of the data points. Every dot represents 200 
data points. 



factor and size of the analyzed features. All these factors con- 
tribute to the scatter between the existing contrast measure- 
ments. 

Our results differ from earlier observations of the contrast 
of bright features, specially when considering magnetic sig- 
nals B/ n > 200 G at disk center. Previous measurements 
of disk center facular contrast have frequently yielded posi- 
tive values, although they usually were close to zero. Thus, 
from multi-color photometric images, Lawrence (1988) mea- 
sures Cf ac ~ 0.005 at disk center, and Lawrence et al. ( 1988[ ) 
find Cf ac = 0.007 ± 0.00 1. In f a ct, ou r results agree better 
with those of Topka et al. ( 1992 1997) and Lawrence et al. 



( |1993[ ) despite the difference in spatial resolution and stud- 
ied wavelength. For 200 G < (B/fi) < 600 G the agreement 
is also surprisingly good. Nevertheless, these authors distin- 
guish between active regions and quiet Sun. For active regions 
they always measure a negative contrast around disk center 
(for /i = 0.97 and [i — 0.99) irrespective of the magne- 
togram signal, while we get slighly positive contrast values for 
B/[i < 200 G in agreement with their results for the network. 
Since at these field strengths most of our signal originates in 
the network, this agreement is probably not surprising. 

Chapman & Klabunde (1982 ) claim that the contrast shows 
a sharp increase near the limb (and even fit a /i -1 dependence). 
We find that Cf ac peaks between /i = 0.5 and fi = 0.2, depend- 
ing on the magnetic strength of the signal, and then decreases 
towards the limb (Fig. g). Libbrecht & Kuhn (|1984 |1985|) also 



find this behaviour; however, they give fj, < 0.2 for the peak of 
the contrast. Wang & Zirin (1987) and Spruit's hot wall model 
also give a similar value for the /i at which Cf ac peaks. It is 
worth noting that Libbrecht & Kuhn ( |1984| , |1985| ) and Wang & 
Zirin (1987) do not take into account the magnetic field of the 
observed feature, which makes the comparison between our re- 
sults and theirs more difficult, as the CLV obtained is different 
when features are selected according to their brightness rather 
than the magnetogram signal. In the former case there is a bias 
towards brighter features. Our results indicate that the higher 
the magnetic signal, the smaller the /i-value at which the con- 
trast peaks (see Fig. |^a) so that network-like features dominate 
at disk center and features with increasingly large Bj [l closer 
to the limb. This should move the peak of the contrast to smaller 
// when the brightest features are searched for, than when mag- 
netograms are used to identify faculae. Finally, it should be 
pointed out that, for increasingly smaller B //j, values the con- 
trast becomes increasingly independent of n; this agrees with 
the conclusion of Ermolli et al. (1999 ) that the network contrast 
is almost independent of /i. 

It is remarkable that an expression as given by Eq. (Q) re- 
produces the dependence of the contrast of bright features on 
their position (/i) and on the magnetic flux per pixel {B/fi), 
within the range 0.1 < fi < 1 and 17G < (B/fi) < 630G. A 
relative accuracy of better than 10% is achieved almost every- 
where within this domain. However, this multivariate analy- 
sis is only a first step and considerable further work needs to 
be done, since two other relevant parameters for the contrast, 
namely the wavelength and the spatial resolution, are kept at 
fixed values (those prescribed by MDI) in our analysis. A 4- 
dimensional data set is thus needed. A first step was taken by 



8 



A. Ortiz et al.: On the contrast of solar photospheric faculae and network elements 




E 2 

o 









(b) - 


;++ +++ 






(x10- 2 ) _ 



a. 

m 



II 2 

o 




200 400 
Magnetogram signal / /j. 



x10' 



600 



Fig. 8. Dependence on the magnetic flux per pixel, B/fi, of: a) 
Mmax; b) C£ c ax times 10 2 ; c) Q° c ax /(B//i) times 10 4 . See text 
for details. 



Lawrence et al. (1993), who compared observations from dif- 
ferent instruments. At least some further progress in this direc- 
tion can be achieved by employing MDI high resolution data, 
although off-center pointing is required. 

4.2. Comparison with flux-tube models 

MHD models including self-consistent energy transfer predict 
that small flux tubes (diameters smaller than 300 km) appear 
bright at disk center but with decreasing contrast near the limb; 
somewhat larger tubes are predicted to appear dark at disk cen- 
ter but bright near the limb, and finally, very large flux tubes 
(pores and sunspots; not considered in this study) are predicted 
to be dark everywhere (e.g., Knolker & Schiissler 1988| , 1989). 
In such models the contrast at /i w 0.1 is largely determined 
by the brightness of the bottom of the flux tube (and the bright- 
ness of its surroundings, e.g. granular down flow lanes), while 
the CLV of the contrast is strongly influenced by the visibility 
of the hot walls. The bottom of a flux tube is defined as the 
horizontal optical depth unity surface in the interior of a flux 
tube. 



Our results are qualitatively in accordance with this predic- 
tion if we make two reasonable assumptions. First, the network 
and facular features are composed of a mixture of spatially un- 
resolved flux tubes of different sizes. Second, the average size 
of the flux tubes increases with increasing magnetogram sig- 
nal or filling factor. Under these assumptions the upper panels 
of Fig. [| refer to, on average, small flux tubes which domi- 
nate the network, while the lower panels of that figure refer 
to larger tubes mostly present in AR faculae. In our study the 
contrast always has a minimum at /i = 1 and increases with 
decreasing /i (as part of the hot wall becomes visible), until a 
maximum when the contrast peaks (the maximum surface of 
the hot wall is seen). Closer to the limb the contrast decreases 
as less wall surface is exposed. There are, however, clear differ- 
ences between small Bj /i network flux tubes and tubes found 
in AR faculae, i.e. regions with large B j '/i. Network tubes are 
bright everywhere on the solar disk and exhibit a low contrast 
(Fig. [|b), but a high specific contrast (Fig. This implies 
that network flux tubes are brighter than AR flux tubes and 
partly reflects the fact that network flux tubes are hotter than 
AR tubes (e.g., Solanki & Brigljevic |1992[ Solanki |1993[ ). The 
greater brightness at large /i implies that network flux tubes 
have a hotter bottom than larger flux tubes. Since this is also 
true at /i < 0.6, it suggests that the walls of smaller tubes, or 
of tubes in regions with lower filling factor, are hotter as well. 
This is in agreement with the theoretical finding of Deinzer et 
al. (1984b) that the inflow of radiation into the tube leads to a 



cooling of the surroundings and a lowering of the temperature 
of the walls. This temperature reduction is indeed predicted to 
be greater for larger flux tubes (Knolker & Schiissler 1988 ). 



A mixture of flux-tube sizes at a given B / fj, is needed be- 
cause the CLV of Cfac at small B//J, does not agree with the 
predictions for any size of flux tube. The model flux tubes are 
all bright over only a relatively small range of /i values. Hence 
the mixture of flux tube sizes is needed in order to produce 
a relatively /i-independent contrast, as exhibited by magnetic 
features at small B//J,. As can be seen in Fig. ||, the contrast 
shows a more pronounced CLV as tube size increases, in accor- 
dance with the hot wall model, and larger tubes have a negative 
contrast at disk center, as predicted. The high specific contrast 
of small B/fi features (Fig. ||c), and the fact that their contrast 
is positive over the whole solar disk indicates that a change in 
the magnetic flux of the network has a much larger contribution 
to the change of the irradiance than a similar change in flux in 
active regions. 

From Fig. ^|a we can determine the heliocentric angles that 
make the contrast peak, max . For the intervals displayed on 
Fig. I 6 max is 63°, 55°, 58°, 62°, 66°, 68°, 72° and 77°, 
respectively. Assuming the hot wall model with a simplified 
cylindrical geometry for the flux tubes, a Wilson depression 
Zyj of 150 km (Spruit 1976) and the derived 6 max values, it is 
possible to roughly estimate the average value of the tube di- 
ameter for each magnetic range. Taking into account that the 
maximum depth of the wall Z\y is seen when the angle be- 
tween the local vertical to the tube and the line of sight is 
equal to the heliocentric angle, then the diameter D should 
be D — Zw tan(# max ). Applying this approximation to our 
observations, we obtain diameters of 290, 210, 240, 281, 334, 
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365, 460, and 650 km respectively, for the mentioned 6* max val- 
ues and their respective magnetic ranges. These diameters are 
estimated to be uncertain by approximately a factor of two. For 
example, uncertainties in Zyv translate into proportionate rela- 
tive uncertainties in D. 

Finally, we wish to draw attention to the wiggle of the mea- 
sured contrasts around fi = 0.95 in Figs. [}]and ^| This can be 
observed at all magnetic strengths. At fi = 1 the contrast has a 
minimum value, then increases, descending a bit later for still 
smaller /i's before finally increasing slowly towards the limb. 
Topka et al. (1992) show some of these variations very close 
to disk center in Fig. 5 of their paper. They argue that such 
variations are partly due to the inclination of the flux tubes of 
opposite polarities toward each other in the active region they 
observe. However, we average over many network elements on 
multiple days and the persistence of such a structure is surpris- 
ing, in particular also for small B/fi values where the statistics 
are extremely good. 

5. Conclusions 

The results presented in this work indicate that the CLV of 
the continuum contrast of magnetic features changes gradually 
with magnetogram signal (or magnetic filling factor), such that 
the contrasts of AR faculae and the network exhibit a very dif- 
ferent CLV, in general agreement with the results of Topka et 
al. (1992, 1997) and Lawrence et al. (1993l. A possible reason 
for the difference is that the populations of magnetic flux tubes 
found in the two kinds of features are, on average, different in 
size, in agreement with the conclusions of earlier investigations 
(Grossmann-Doerth et al. 11994 Keller 1 19921). 



Stronger magnetogram signals, corresp onding to wider flux 
tubes on average (Grossmann-Doerth et al. 1994 ), appear dark 
at disk centre, but bright near the limb, while the weakest 
signals (on average narrower flux tubes) are almost equally 
bright at disk centre and near the limb. This result is in good 
agreement with t he predictions of theoretical flux tu be mod- 
els (Deinzer et al |l984a| |l984b[ Knolker et al. |1988| ; Knolker 
& Schiissler 1988) if there is a distribution of flux tube sizes 



present within an MDI pixel. Because network elements are 
bright over the whole solar disk their contribution to irradi- 
ance variations is significant and needs to be taken into account 
when reconstructing variations of the total solar irradiance. 

One advantage of the present investigation relative to that of 
Topka et al. ( |1992[ |1997[ ) is that by using full disk MDI data we 
have a result for a very well defined spatial resolution, so that 
any models derived on the basis of these results can be directly 
used for reconstructing total and spectral solar irradiance, as for 
instance measured by VIRGO (Frohlich et al. 1995), without 
further adjustment. 

A new result of this work is that, with a simple expression, 
we can predict the contrast of a bright feature, from network 
and small tubes to faculae of different sizes, given its position 
and magnetic strength within a certain range, and reproduce 
simultaneously the Cf ac (/i) and Cf ac (£?//i) dependences. 

In a next step the dependence of the contrast on wavelength 
(for given ji and B //i) must be determined, as well as the de- 
pendence on spatial resolution. The later dependence is of par- 



ticular interest also because the investigations of Lawrence et 
al. (1993) and Topka et al. (1997) give similar values of the 
contrast near /i = 1 as we find, although the spatial resolution 
of the La Palma data employed by these authors is almost an 
order of magnitude higher than that of the MDI full disk data 
(0.5"versus 4"). Closer to the limb Topka et al. ( |1997[ ) obtain 
contrasts a factor of two higher. Whether this is due to the dif- 
ferent wavelengths observed or has another source needs to be 
investigated. 
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Fig. 3. Facular and network contrast as a function of n for eight intervals of the strength of the magnetic field, from network 
values (top left panel) to strong faculae (lower right). A dashed line at Cf ac = has been plotted. The solid curves represent a 
second order polynomial least-squares fit to the points. Every dot represents 40 data points, /j, = 1 is the disk center; \i = is the 
limb. 



12 



A. Ortiz et al.: On the contrast of solar photospheric faculae and network elements 



0.08 



m 0.04 
o 



o 
o 



0.00 



0.08 



m 0.04 

o 




o 
o 



0.00 



0.08 



w 0.04 

o 



o 
o 



0.00 



0.08 



m 0.04 

a 



o 

o 



0.00 



0.92</x< 1 .0 



0.78</x<0.82 



0.54</z<0.66 



0.32</x<0.4£ 




0.88</x<0.92 



0.66</x<0.74 




0.42</uK0.58 




0.20</x<0.40 




200 400 600 

Magnetogram signal//^ 



800 200 400 600 800 

Magnetogram signal///, 



Fig. 4. Dependence of the contrast on the absolute value of the magnetogram signal, corrected for foreshortening effects. The 
solar disk has been divided into eight bins, from center to limb. Note that some ^t-bins overlap (see text for details). As in Fig. |l 
a dashed line at Cf ac = has been plotted and solid curves represent a second order polynomial regression. Every dot represents 
40 data points. 



